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Gamma-ray Burst UV/optical afterglow polarimetry as a 
probe of Quantum Gravity 
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ABSTRACT 

A possible birefringence effect that arises in quantum gravity leads to a frequency- 
dependent rotation of the polarization angle of linearly polarized emission from distant 
sources. Here we use the UV/optical polarization data of the afterglows of GRB 020813 
and GRB 021004 to constrain this effect. We find an upper limit on the Gambini & 
Pulin birefringence parameter \tj\ < 2 x 10~ 7 . This limit is of 3 orders better than 
the previous limits from observations of AGNs and of the Crab pulsar. Much stronger 
limits may be obtained by the future observation of polarization of the prompt 7-rays. 
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Modifications to the standard dispersion relation of high 
energy particles arise in various theories of quantum gravity 
(Amelino-Camelia et al. 1998; Amelino-Camelia 2000; Piran 
2005). These modifications result generically in a modifica- 
tion to the standard dispersion relation. This, in turn, in- 
troduces an energy dependent photon's velocity. Attempts 
to explore this possibility and to set limit or observe it 
are typically based on time of flight tests of very high en- 
ergy photons. In specific cases and in particular in loop 
quantum gravity, the deviations from the standard disper- 
sion relation take a particular form in which photons with 
right handed and left handed polarization travel with dif- 
ferent velocities (Gambini & Pulin 1999). This leads to a 
frequency dependent rotation of the polarization vector of 
a linearly polarized light. Observations of a linear polar- 
ization can, therefore, set limits on the magnitude of the 
quantum gravity birefringence parameter (Gambini & Pulin 
1999; Gleiser & Kozameh 2001; Mitrofanov 2003; Jacobson 
et al. 2004; Kahniashvili, Gogoberidze & Ratra 2006). The 
presence of linear polarization in the UV/optical spectrum 
of distant galaxies and in the the X-ray spectrum of the Crab 
Nebula sets upper bounds on the birefringence parameter: 
77 < 10~ 4 (Gleiser & Kozameh 2001; Kaaret 2004). Much 
smaller limit, r/ < 10 -14 , can be obtained using the linear 
polarization detection in GRB 021206 (Mitrofanov 2003; Ja- 
cobson et al. 2004). However, the original claims of a detec- 
tion of polarization of this burst (Coburn & Boggs 2003) 
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have been refuted convincingly (Rutledge & Fox 2004; Wig- 
ger et al. 2004). The polarization analysis of GRB 030519B 
also yielded null result (Wigger et al. 2004). Possible ev- 
idence for linearly polarized 7— rays in GRB 930131 and 
GRB 960924 (Willis et al. 2005) was also reported. If cor- 
rect it will provide a comparable upper limit. However, it 
is not clear yet how robust are the novel methods used to 
estimate the polarization in these bursts. We explore here 
the implications and limits that can be set by observations of 
linear polarization from the optical afterglow of GRB 020813 
(Barth et al. 2003) and GRB 021004 (Lazzati et al. 2003). 

Gambini and Pullin (1999) have shown that within 
loop quantum gravity photons with left and right circu- 
lar polarizations will travel at different velocities v± = 
c(l ± rjhcu/E pl ), where ± denote the different polarization 
states and r\ is a dimensionless constant. Here we discuss 
the generalized form (see also Mitrofanov 2003 and the ref- 
erences therein) 



v± = c[l ± r](hLu/E pl 



(1) 



where r\ 1 ' n characterizes the scale of breakdown of the 
standard dispersion relation in units of the Planck's energy 
E p] = 1.22 x 10 19 GeV. Naturally, one would expect 77 ~ 1 
corresponding to a typical energy scale of the Planck energy 
or slightly higher, while a very low limit on r\ would render 
such a modification unlikely. 

Since the linear polarization is a superposition of two 
monochromatic waves with opposite circular polarizations, 
the plane of linear polarization is subject to a rotation along 
the photons' path because of the difference between the 
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Figure 1. F(z,n) in the rcdshift range 0.01 — 10 for Cm = 0.28, 
n A = 0.72, and h = 0.73. 



two circular components (Gambini & Pulin 1999; Gleiser 
& Kozameh 2001; Mitrofanov 2003; Jacobson et al. 2004): 

^ = ^)«+i^W, (2) 

C lr> 



where l v = ^JfiG/c 3 — hc/E pl 
Planck's length scale and 

dz 



1.6 x 10' 



cm is the 



hHo (1 + 2) v /O a + «m(1 + z) 3 



is the differential distance the photons has travelled (for a 
flat cosmological model, i.e., 12m + = 1), h ~ 0.73 is the 
current Hubble constant in units of Ho = 100 km/sec/Mpc. 

The observed polarization angle at frequency v ^ s of a 
photon emitted at a redshift z with an intrinsic polarization 
angle $0 is 



$ (n) = $0 + 7.94 x 10 60 v ^ TT (2nu obs ) n+1 F(z,n) 



(3) 



The function F(z, n) depends on the cosmology and can be 
expressed as 



F(z,n) = ( — ) 
v ' ; ^0.73 ; 



(l + z) n dz 



v /«a + «m(1 + ^) 3 



(4) 



The function F ~ 1 for 2 = 1. F(£, 1) a z for 2 <C 1. 
F(z, 1) tx 2 1/2 for z > 1, while F(«,2) <x 2 for 2 < 1 and 
F(z,2) oc 2 3/2 for 2 > 1. The slow increase of F(z,l) for 
2 > 1 suggests that there is no advantage in going to very 
large z values. 

For n — 1, we have A<6m — $(i) — $0 = 
0-9rj( 1 q°2 I | Iz ) 2 F(z, 1). Therefore, the relevant frequencies for 
observing the birefringence effect are u b s > v^ ir = 1.05 x 
10 12 7?~ 1//2 F(2, 1)~ 1/2 , where vu r is defined as A$ (n j (vur) = 



1. For n = 2, A$ 



(2) 
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$0 = 0.3r/(- 



5 F(2,2). So 



> 1022 Hz > 

we have to probe the birefringence effect with polarimetry 
of the ~ 50 MeV photons from a source at 2 ~ 1. This is 
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quite difficult because usually at such high energy, there are 
not enough photons to allow a reliable polarization measure- 
ment. Consequently, there was no X— ray/7— ray polarime- 
ter planed for energy higher than 30 MeV (McConnell et al. 
2006). In this work we focus, therefore, on the case of n = 1. 

At u b s <C fbir there is no significant observable effect. 
On the other hand at u a b a ^> Vbir the rotation is so large 
that it will erase any intrinsic linear polarization (Gleiser & 
Kozameh 2001; Jacobson et al. 2004). Therefore an detection 
of high linear polarization at a frequency v \, s from a source 
at 2 puts an upper limit on r\: 



1 V<( 



Vobs 
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) F(z, 1) , or?? < 3( 
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)- 3 F(2,2)-\(5) 



In a limited frequency range around Vbir we expect to ob- 
serve the birefringence effect as a frequency dependent lin- 
ear polarization vector. Such an observation could confirm 
the existence of a birefringence effect and the corresponding 
quantum gravity corrections provided that we can rule out 
an intrinsic origin for the rotation of the polarization vector. 

The UV/optical spectropolarimetry of GRB afterglows 
has been reported in several events (Covino et al. 1999; Wi- 
jers et al. 1999; Barth et al. 2003; Bersier et al. 2003; Greiner 
et al. 2003; Rol et al. 2003; Gorosabel et al. 2004). Good 
quality results were obtained on GRB 020813 (Barth et al. 
2003) and GRB 021004 (Lazzati et al. 2003). We use these 
two bursts to constrain the possible birefringence effect. 

GRB 020813 was detected by High Energy Transient 
Explore 2 (HETE-2) on August 13rd 2002 (Villasenor et 
al. 2002) and found at the redshift 2 ~ 1.255 (Barth et al. 
2003). The UV/optical [(4000, 5000, 6300, 7300, 8300) ± 
500A] polarization measurements of the afterglow of GRB 
020813 were carried out at ~ 5.16, 6.27, and 7.36 hours 
after the burst, respectively. The observed polarization is in 
the range of 1.8% - 2.4% (see Table 2 of (Barth et al. 2003) 
for details) and the scatter is very small, in particular for 
the data collected at a given time. At t ~ 5.16 hours, the 
corresponding polarization angles are (161 ±1, 159±1, 158± 
1, 155 ± 1, 155 ± 1) deg, respectively. At t ~ 6.27 hours, the 
corresponding polarization angles are (160±1, 155±1, 151± 
1, 150 ± 1, 151 ± 2) deg, respectively. At t ~ 7.36 hours, the 
corresponding polarization angles are (153±1, 149±1, 156± 
1, 153 ± 1, 149 ± 1) deg, respectively (Barth et al. 2003). 

To constrain 77 we look for a similar frequency depen- 
dent trend in this data. We allow a variation in the po- 
larization from one time to another and consider only the 
relative polarization angle. To do so we subtract, at each 
time, the average polarization angle (which can vary intrin- 
sically with time) and consider only the polarization angle 
with respect to this average (see Fig. 2). The scatter in the 
shift between different times provides an estimate of the er- 
ror in the data. The best fitted r\ — 1.9 x 10~ 7 and r\ — 
is marginally consistent at the 2a level (see Fig.2)t Upper 
limits at the 3<r level that can be drawn from this figure are: 
-6 x 10~ 8 < 77 < 1.2 x 10" 6 . These limits are two orders 
of magnitude better than the limit obtained using optical 
polarization data of one distant galaxy (Gleiser & Kozameh 



f The direct time-averaged polarization data yields a less tight 
constraint on the birefringence effect. In such a treatment, 77 = 
is consistent with the data at lc level. 
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Figure 2. Fit to the afterglow polarimetry data of GRB 020813 
(the upper panels; the angles are time-averaged) and GRB 021004 
(the lower panels). The data are taken from (Barth et al. 2003) 
and (Lazzati et al. 2003), respectively. Left: Linear fit of the po- 
larization angles vs. the quantities of v^ bs /W 2 ®, taking double 
errbars into account. Right: Joint confidence intervals (la, 2a, 
and 3a) in the <I>o — J? plane. The cross indicates the best fit, 
corresponding to the minimal chi-square value X 2 i n - 



2001) and X-ray polarimetry of the Crab Nebula (Kaaret 
2004). 

GRB 021004 was detected by HETE-2 (Shirasaki, 
Mathewson & Ford 2002) on October 4th 2002 and pinned 
down at the redshift z ~ 2.328 by follow-up observations 
(Mirabal et al. 2002). Here we directly take the reduced 
spectropolarimetry data presented in Lazzati et al. (2003) 
using the Very Large Telescope of the European Southern 
Observatory. As shown in the bottom left panel of Figure 2, 
we find the polarization angle of this burst is proportional to 
^obs with a negative rather than a positive coefficient. The 
best fit indicates r] — —1.1 x 10 -7 . The dataset is consistent 
with the possibility of no shift at all (i.e., r\ — 0) at the 2a 
confidence level. 

Now we combine GRB020813 and GRB021004 to fur- 
ther investigate the possible birefringence effect that arises 
in quantum gravity, that is, mainly to constrain the param- 
eter r/. To realize this, we apply the Chi-Square fitting to the 
data of both bursts and marginalize the Chi-Square merit 
function over all the possible polarization angles in order to 
derive the probability only for the parameter r\. The pro- 
cedure is to (1) calculate the probability for each given set 
r/ and polarization angle for either burst; (2) integrate the 
probabilities over polarization angle and then derive two av- 
erage probabilities for each r\ for either burst; and (3) derive 
the joint probability for each r\ for two bursts by using the 
formula Pfo) oc Exp[-(xo208i3,» 7i + Xo208i3,^ )/% where i 
means the i-th r?. GRB020813 and GRB021004 are treated 
as independent events during the joint analysis. Shown in 
Figure 3 are the likelihood functions for the parameter 77 



Figure 3. Individual and joint constraints on rj from the 
UV/optical polarimetry of GRB 020813 and GRB 021004. For 
the joint case, —2 X 10 — 7 < rj < 1.4 X 10~ 7 at the 3cr confidence 
level. 



of GRB020318 (blue dash curve), GRB021004 (olive dash 
curve), and their combination (red solid curve). At the 3<r 
confidence level, we have -2 x 10~ 7 < r\ < 1.4 x 10 -7 . This 
limit is about 3 orders better than the previous estimates. 
For such stringent constraint on rj, the birefringence effect 
on our initial polarization measurement in a limited band- 
width (500A for GRB 020813 and 200A for GRB 021004) is 
unimportant and can be ignored. 

It should be noted that the rotation of the polariza- 
tion plane can also be induced by Faraday rotation, so it 
is necessary to test whether the Faraday rotation could af- 
fect the results discussed in this paper. It is well known 
that the rotation angle induced by the Faraday rotation is 
Act = 8.1 x 10 5 A^ J n e B\\dL, where A m is the wavelength 
in units of meter, n e is the electron number density per 
cm 3 , Bu is the magnetic field strength (Gauss) parallel the 
propagating direction, and L is the distance in units of pc. 
Typically, with A m ~ 10~ 6 , n e ~ 10" 4 , By < 10~ 12 G and 
L ~ 10 11 , we have Ao < 10~ n . As can be seen, even for 
z ~ 10, the rotation angle at UV/optical and higher band is 
very small, therefore the Faraday rotation is insignificant. 

A detection of polarization of the prompt emission of 
GRBs (at such high frequency) will set immediately a much 
stronger limit on 77 (Gambini & Pulin 1999). As mentioned 
before, a possible evidence for the linearly polarized 7— rays 
(~ a few hundred keV) has been reported (Coburn & Boggs 

2003) and ruled out (Rutledge & Fox 2004; Wigger et al. 

2004) in GRB 021206. Wigger et al. (Wigger et al. 2004) 
also suggest that there is no evidence for polarization of the 
prompt 7-ray emission in GRB 030519B. On the other hand 
polarization was detected in GRB 930131 and GRB 960924 
(Willis et al. 2005) and possibly in GRB 041219a (McGlynn 
et al. 2006). These results are consistent with each other as 
(apart from the refuted claims of (Coburn & Boggs 2003)) 
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the errors are very large. The situation is inconclusive and 
additional data is needed before the corresponding strong 
limits on rj that have been suggested on the basis of this 
burst (Jacobson et al. 2004) can be accepted. The planned 
X-ray/7— ray polarimitors POLAR (Produit et al. 2005) or 
XPOL (Costa et al. 2006) could resolve this issue in the fu- 
ture. Even such polarimeters won't help with a significant 
limit on n = 2. To reach the limit 77 = 1 for n — 2 one needs 
a polarization measure of 10 — 100 McV signal from a cosmo- 
logical distance. However, even for the bursts with a hard, 
highly-polarized component as energetic as the one detected 
in GRB 941017 (Gonzalez et al. 2003), such a polarimetry 
requires a very large area (~ 1 m 2 ) which is beyond the 
goals of current design. 
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